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Forced flow of a Fermi liquid is studied for a cell geometry consisting 
of two planes with a separation on the order of the mean free path. An 

approximate transport equation is used to derive an integral equation for the 

velocity profile, which is solved numerically. Results for the total flux 
through the cell, t/hich determines the dissipation, are given as a function 

of the Knudsen number x, (ratio of cell thickness to mean free path.) Effects 

of specular reflection af the boundaries are considered. It is found that the 
dissipation has a minimum at x. 2! H, and behaves linearly for x. * 3. Implica- 
tions for present experimentation are discusser. 
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ABSTRACT 

Forced flow of a Fermi liquid is studied for a cell geometry consisting 
of two planes with a separation on the order of the mean free path. An 

approximate transport equation is used to derive an integral equation for the 

velocity profile, which is solved numerically. Results for the total flux 
through the cell, which determines the dissipation, are given as a function 

of the Knudsen number v. (ratio of cell thickness to mean free path.) Effects 

of specular reflection at the boundaries are considered. It is found that the 
dissipation has a minimum at — 35 , and behaves linearly for H > 3. Implica- 
tions for present experimentation are discussed. 
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I. Introduction 

Recent experiments 'fove provided precise measurements of the transport 
properties of the various phases of ^He to below 1 mX a . In this temperature 
regime the mean free path of the quasiparticles approaches or exceeds the 
dimensions of the experimental apparatus. This means that quasiparticles 
scattered from one boundary of the apparatus will not relax into local equili- 
brium with the rest of the fluid before they reach another boundary. The mean 
free path in the normal phase may approach 50 ja, and is much larger in the 

4 

superfluid phases , while Andronikashvili cells presently in use have a thick- 
ness on the order of 100 jm. In consequence, ch* viscosity can no longer be 
related to measured quantities (e.g. drag forces) by the usual formulas of 
macroscopic fluid or superfluid mechanics. It is therefore necessary to under- 
stand such finite size effects in order to interpret present experiments. 

This paper aims to shed some light on this problem by analyzing a parti- 
cular case, namely Poisleulle flow, or forced flow of a normal fermi liquid 
between stationary infinite planes. We will calculate the total flux through 
the cell, which determines the rate of energy dissipation. In the large-cell 
limit this flux is inversely proportional to the viscosity, so the finite-cell 
flux may be taken to be irversely proportional to an "effective viscosity" that 
depends on the cell thickness. The ratio of this "effective viscosity" to the 
true viscosity of the bulk fluid conveniently characterizes the effect from an 

intuitive standpoint, although the dissipation is more closely related to 

1 3 

the measured quantities. Some recent experiments ’ will be considered in 
the context of this theory, and possible relevance to future experiments 
discussed in the conclusion. 

II. Choice of a Model 

Consider a Fermi liquid flowing between stationary walls. The flow is 
sustained by a static, uniform force and is taken to be in the x-direction. 
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The cell is taken to be of infinite extent in the x- and y-directions and to 
extend from z ■ -D/2 to z ■ D/2. (See fig. 1) The system clearly has mirror 
symmetry in the xy-plane. 

We will consider the quasiparticle transport equation for this system. 

The Landau Interaction between quasiparticles^ will be seen to cancel out of 
this equation, and enters the prob’em only through the quasiparticle lifetime 
T^. We will treat the latter as a given quantity, since we are Interested 
in the comparison of a finite and an Infinite sample of the same fluid. If 
we treat. particle-particle scattering in the same manner in both cases, the 
deficiencies of the relaxation time approximation should not drastically affect 
our basic results. We will see that r contains the entire temperature as well 
as interaction dependence in the problem. 

The simplest model which can be used for the scattering of particles from 
a wall is that of completely diffuse reflection. This means that the particles 
emerging from the wall after colliding with it possess a distribution which 
is isotropic in the reference frame of the wall and which corresponds to local 
equilibrium conditions at the point of collision. A more general model, which 
we shall use, is to assume that a fraction e of the particles are specularly 
reflected^, that is, the component of their momentum perpendicular to the wall 
is reversed and the other components are unchanged by scattering from the wall. 
This is certainly not the most general possible reflection process, but since 
very little is known about the quasiparticle-reflecting properties of experimen- 
tal surfaces, it should be sufficient for our purposes. We consider Z as an 
edjustable parameter; we will see that the 2 -dependence of our results can 
easily be separated from the other dependences. 

III. Derivation of the Integral Equation 

The number of quasiparcicles per unit phase volume with momentum p at a 
height z in the cell is written n~(z). In the absence of the applied force 
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field this becomes equal 10 Che global equilibrium distribution function, 
which is written 



e^P 


■U> + 


( 1 ) 


. 0 2 ,* * 

Here 3 ■ i/k_T, € ■ p /2m whers m is the effective mass, and ^ is the 

8 p 

chemical potential. We make the approximation chat Che temperature and 

g 

chemical potential remain constant and uniform . The average velocity of 
the fluid at height z, or "velocity profile", is given by 


u(z) “ 5 v- n-(z) , 

p P P 

— e — * 

where v^ is the velocity of a quasiparticle with momentum p. Since Che fluid 
flow ij entirely in the x-direccion, we may write 


u(z) 


P x 

£ ~ n p< z) 

pm r 


( 2 ) 


we also define 


5n-(z) ■ n-*(z) - n^ , (3) 

P P P 

s p< 2 > * £ p + S, fin p' U) » 

where f^i.is the Landau quasiparticle interaction function. We also define 
the local equilibrium distribution function and the deviation of the actual 


distribution from it: 



s 




a 3(e-(z)-P x u(z)-u) +l 


Sn-(z) • n-(z)-n- (z) . 
P P P 


Wd observe from (2), (5), and (6) chat 


' p6n-*(z) ■ 0 . 

P P 


The transport equation is given by 


t~(z) + 7 r s r 7 r n r<*) + Pr(z) , 7-n- > (z) • i(n 1 ] , 
dc pprp p pp p 


where Itn^i] denotes the collision integral and the force F-(z) on a quasi- 
particle with momentum p includes both the external field F and a mean field 
part : 


F-(z) • F • V- J f~, 5n-, (z) . 

P r pi PP P 


We now assume that the system is in a steady state, and approximate Che 


collision integral by 


I[n-,] “ 


6n-(z) 


We also linearize in the quantities F, u, fin", &n—, and use the fact that F 

P P 

is in the x-direction while all spatial gradients are in the z-direction. 

We find that the term containing f-*-*, which comes from the third term on the 

PP 

left side of the transport equation (6) is cancelled by a similar term arising 
from the second term in (6). The transport equation becomes 
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'_L 





(11) 


This differencial equation in z has Che general solution 

z in° 

*“p<«) - exp (-**z/P z T){a^ : dz , exp(a*z , /p z T)p - r L )(^)} (12) 

r -D/z r z p 

where A- is an arbitrary function of p. When we apply boundary conditions 
to (12) we will obtain an explicit expression for 5n^*(z) as a functional 
of u(z). Requiring that this functional satisfy (7) will then give us the inte- 

A 

gral equation for u(z), from which the flux u u(z)dz is determined. 

We now consider the boundary conditions. If a fraction s of the quasi- 
particles incident at a wall are specularly reflected, then the distribution 
function at the lower wall obeys 


n-(- §)-<!- Ore ♦ |). > 


(13) 


where p* ** (p ,p ,-p ). At the upper wall the distribution function obeys 
x y z 


n-(D - (l - e)«j°+ «-„(§) , p 2 > 0 


(14) 


since n- is independent of the sign of p . To express these conditions in 
P z 

terms of the deviation from local equilibrium, we go back to the definition 
(5) of the local equilibrium distribution function and expand to linear order 
in the perturbation F: 


"p'v- f) - “p + (sjJH5.‘E’ 6 v(' 2 ) 


/ dm 

cV W 


(15) 


Substituting nj(- §) - n^(- f) + §) into (13) and using (15), we find 
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5n-*(z) - efin- (*) 

P P* 




PP 


• 6 v(- ?)} 


(16) 


Consider Che Fermi liquid inceraction terms in (16). Alchough Che'/ 
are formally of che same order as the ocher cerms, our model of vhac happens 
ac che boundaries is sufficienCly crude chat we can probably omit them wlthouC 
making Che description any lees accurate. Recall chac % is a phenomenological 
parameter Chat reflects Che degree of correlation between Incident and 
reflected quasiparticles. Ue assume that che effect of the interaction cerms 
in (16) may be adequately approximated by choosing a different value of Z 
than we would have done had no interactions been present. Even without this 
assumption we have no way of calculating e for the surfaces of interest, and 
retaining the interaction term would make the calculation much more difficult. 
Thus (16) becomes 

4 "p(* 2) ' §) ‘ (a^) (1 - s,1 VV £)■ p s > o • < l7 > 

# 

The analogous equation for the other boundary is 

56^1) - (n*)a-«>p x »(i) • <18) 

p 

To apply these conditions, we go back to the general solution (12) of 

the transport equation and evaluate it a z ■ -%D and z m -b^D, substituting 

these results into (17) and (18) respectively. This gives us two equations 

for the two unknowns A- and A-. . Putting the solution for A- into (12) gives 

P P* P 

us 6n-(z) for p >0: the solution for p < 0 is found by exploiting the 

P 2 2 
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mtrror symmetry of the problem: 


5n-(z) - 5n-(-z) 

P P* 

and 

u(z) - u(-z) . 


When we follow this procedure, the results are 


•:exp(-ra*D/2p t) 

— . ■ 2 


. , -«*Z/D,T| ■‘‘•‘Kf) 

5n-(z) - e r z 

p 'sexp(-m*D/2p T)-exp(m*D/2p T) :exp(-m' r D/2p T)-exp(m*D/2p t) 

z z z z 

D/2 


* Bn 


■D/2 


-D/2 


(21a) 


iiHt) - 

P 


exp(-n*D/2p t) 

— i- z . . 


'-:exp(m*D/2p "1-exp (-ra*D/2p r) -;exp(m*D/2p T)-exp (-m*D/2p r) 

- z z z 

in 0 

; di’«xp(»*.-/p J D(|a,. i)+/d»'«p<«*,Vp i T>(|li,- ftyrJ-JOj, ? < o . 

- z p 

(21b) 


D/2 


-D/2 


-D/2 


It will be noted hat the distribution function is discontinuous at p ■ 0. 

z 

This is because the particles which have just been scattered by a wall have 
a substantially different distribution than those about to be scattered due 
to the large change in momentum that occurs in scattering. 

We can now impose he condition (7) that the average of the momentum 
over the deviation fron local equilibrium vanishes. The only component that 
does not v?.nish trivially is the x-component, thus the equation of interest is 


0 * j d3 P P x 6 ”-(z) 


( 22 ) 
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Inspection of (21 a) And (21b) shows that (dn /de ) is the only factor 

P P 

in 3n^(z) that varies rapidly as a function of |p| near p p . Thus, as long 

as T « T pl we may replace it by -6 ( 2 p “ '-p^* Uh« n w« integrate using spherical 

coordinates in (22), the effect of this 6-function is to set |p| ■ p everywhere. 

r 

The y~ integral is trivial, so that, defining ^ - cosd and the mean free path 
l ■ p p T/m*, (22) bee ?mes 


I* 1 . f . 2. -z/lif “(l-€)u (f) 

- I 4a(1-u )• 

* a L . t — / m . v 




0 'Sexp(-D/24a)-«xp(D/24i) -:exp (-D/2Aa )-exp CD/ 2 4a ) 

D/2 z 

'dz'exp(z , /Ai)(|“i * ^•)+J‘dz , exp(z'/iA)(|j,- f)} 


j 

-D/2 


-D/2 


+ (l*u 2 )exp(-z/Aa)'[- 


-(l-s)u^ 


exp(-D/24u) 


exp(D/2^)-exp v D/2ia) eexp(D/24i )-exp(-D/24a 


D/2 


j a*' exp (z’/fc )(!£,- dz , exp(z , /4 i )(|j,- £■)} . 

-D/2 z " D/2 


(23) 


We next integrate by parts wherever du/bz' appears. It is found after 
some algebra that the surface terms from these integrations cancel the terms 
in (23) containing u(D/2). If we define a dimensionless velocity 


U(z) 


m*u ( z ) 
FT 


and divide by -Ft/m* we have 


1 - , sexp(-z/^i)exp(-D/2Ai) . D/2 

0 » j* 4adTi ){-2U(z)- y-; dz'exp(z74 k )(l+U(z')) 

0 ' €exp(-D/24i) - exp(D/24a) ^-D/2 

D/2 


exp(z/4l)exp(D/24i) . 

— 0 ’dz'exp(z7li)(l+U/z')) 

€exp(-D/2ii)-exp(D/2il) _ D / 2 



exp(z' 4 a) 
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•xp(-t'ii) z 

f 4*' •xpU'/fc) (1-HJ («'))♦ 

4a -b/2 


4 * 


•z 

dz'expU'/^Xl+UCz'))* . 

-D/2 j 

(24) 


Further manipulations permit u* co case this in the form 


U(z) 1 • • *xp(*z/4a)| * dz'exp(z'/ 4 a)(1>“HJ(*’ ))— :«xp(-D/4a) 

b * l-:«xp(-D/4* -D/2 

f dz 'exp (z '/ 4* ) (l"HJ(z' )) -*-ex? (*/ 4a )l dz 'exp (z' / 4a ) (l+U(z ' ))-*;exp (-0/ 4a) 

D - /2 J ^b/2 

-z 

J* dz'exp(z'/ 4 a)(1*HJ(z' )) j} • (25) 

D/2 

9 

We now define Che Kr.udsen number (ratio of call thickness to mean free 

?*ch) 




1 

l 



(26) 


a dimensionless coordinate l • z/D, and redefine the functional dependence 
of the dimensionless velocity U so thae U(z) - U(£). With a few more steps 
we obtain 


l 


,1 


uo - ^(1 +o;4a(- *u) 


i 


L-*:exp (-kAa ) 


dCexpC-HjC-C |Aa)(1+uC')) 


(27) 


This is the integral equation for the velocity profile. 

IV. Solution for Flow Rate 

We have solved this equation numerically by doing a Gaussian integration^ 
on^, using a finite element method to convert the integral equation into a 
matrix equation, and then using an iterative technique. This gives a discretized 
representation for the solution U(') which we integrate numerically across 
the cell. We then have the quantity 


h 

j'u(C)dC - <o> - 

A ^ 


p r <u> 


X‘ 


FD 


(28) 


as a function of h and €. In order to see the X-dependence of (u) more 
clearly, we divide (28) by x and plot the result for several values of z in fig. 
2 . 

.ae resulting curves have several striking features. For x > 3 the various 

values of e give a family of straight lines whose slope with respect to x Is 

close to the ratio S/12. We shall see that this is the exact value required 

to recover classical bulk fluid behavior in the limit of large x: one can 

f O 

also verify analytically chat the solution of the Integral equation (27) has 

just this behavior. Thus, the only manifestation of a finite size effect for 

x * 3 is the x-independent offset of these lines. 11 

At small x the quantity (U)/x diverges logarithmically. This divergence 
arises from the rather artificial assumption that quasiparticles scattered 
almost parallel to the wall move slowly away from it and do not Interact with 
it any more. These particles will cake longer than most others to cross the 
cell, and hence will be acted upon by the driving field for a long time, giving 
a large contribution to the total current; however, these particles are the ones 
most likely to scatter against ocher particles before reaching the ocher side. 

A logarithmic divergence of (U)/x with x. can be shown to result from this 
situation. In a real cell, in which quasiparticles reflected parallel to a 
wall will be re-scattered by surface Inhomogeneities, we expect a large but 
finite flux at x * 0. Suppose the inhomogeneities in the surface have typically 
an amplitude a and a wavelength L. Then Che minimum possible angle between 
the wall and the momentum of a quasiparticle traveling away from the wall is 
approximately a/L. The quantity (U)/x will cease going as - M x when x ^ a/L; 
instead, we expect it to go to a constant value of roughly *i(L/a) . 

In the intermediate x region we find that (U)/x goes through a minimum 
in the vicinity of K ■ i, and then approaches the linear behavior seen at large 
x. If we extrapolate the high-x behavior down into the intermediate region, 


and find the difference between this linear function and the actual (U)/h, w« 
obtain a measure of the deviation from linear behavior. This deviation, which 
we call 6, la displayed in Table 1; it la seen to drop off rapidly as H increases. 

The e-dependence of the results is also worthy of note. It can be well 

approximated by an additive term Independent of h. This implies that mr 

choice of e will not affect the qualitative structure of tfc* results, 'we note 

that a very similar result has been found for the case of classical rarlfied 

. l 2 

gas dynamics 

The flow rate (u) through the cell determines the rate at which the applied 
force F does work on the system; this dissipation is what will usually be 
measured in an experiment. One way to characterize the finite size effect is 
to say that the fluid behaves as if it had a different viscosity than would 
an infinite sample of the same fluid. We define this "effective viscosity" 
to be that which would give the correct flux (or dissipation) of the actual 
cell when substituted into the equations of hydrodynamics. In the classical 
theory the flux is 


o 


<U> 


D^Fn 

12r| 


O 


(29) 


where n is the number density in the fluid. For the finite cell we define 



(30) 


We now divide both sides by the viscosity of the bulk fluid, using (28) and 

o o 


the relation 

C 


•• Vlk " 5 p fV t * 

and obtain 

^eff J_ * 2 

Vu 12 w 


O 

o 


o 


o 


(31) 


o 


(32) 
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This ratio must obviously approach unity at large K : the asymptotic behavior 

mentioned above for the ratio (u)/k clearly Insures that thi;. will be the case. 
The quantity (32) is plotted in fig. 3 i a function of k for several values 
of •: . 

V. Application to Experiments 

These results have been applied to the interpretation of some experiments 
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performed by Parpla ’ and Archie . In these experiments, a disc-shaped 

3 

cavity in an epoxy cell is filled with He and caused to vibrate torsionally 
about its central axis. The drive voltage which is required to maintain 
constant amplitude with changing temperature is then measured. The measured 
voltage is proportional to the total rate of energy dissipation in the cell, 
which is the sum of the dissipation in the liquid and a "nuisance damping" 
associated only with the apparatus. 

Some results of these experiments are given in fig. 1 of reference 1. The 

effective viscosity, obtained from the driving voltage required to maintain con- 

2 

stant amplitude, is plotted as a function of T ; it is linear above about T ■ 

3.5 mK , and shows an upward curvature from this temperature down to the transi- 
tion temperature (below which it jumps sharply upward.) The present calculation 
was originally undertaken in the hope of understanding this curvature. 

To apply the results of the previous section of this paper, we make 
several observations: (i) the frequency in these experiments is less than 

1 KHz, permitting a steady-state analysis, (ii) the cavity radius is about 
forty times its thickness, permitting an element of the cell to be treated 
as an element of an infinite slab executing translational vibrations with an 
amplitude of r0 Q , where r is the distance of the actual element under considera- 
tion from the axis of the cell, and 9 q is the amplitude of the tcrsional vibra- 
tion, (iii) Edge effects may also be neglected, (iv) the vibrational amplitude 
is very small in these experiments (less than 10 ^ radians) so that a ' -nearized 
theory is well justified. 
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Ic Is easy to relate our parameter K to the < jmperature. We have * ■ D/X ■ 
D/VpT, and In the normal phase the relaxation time is Inversely proportional 
to the square of the temperature.^ Writing 


T - AT* 2 , (33) 

we have 

* . fiii 

v p A • <*♦> 

Since v_ and A depend only on the pressure, we see that the Kr.udsen number is 

f 

proportional to the square o £ the temperature for constant D. The constant 
of proportlorality is found from experimental data; see Table 2. 

We must now relate the dissipation to the dimensionless flow rate found 
in the previous section. In the frame of reference of the wall, there is 
a fictitious force of -ma per particle, where a is the Instantaneous accelera- 
tion of the v;all in the laborscory frame. We use our static theory to obtain 
the response to this force; the relative error involved in doing this will 

-3 

be of order j;t 3s 10 in our case. The dissipation in an element dA of the 
cell is 


dE 


E ■ -nma(u) DcL 1 “ — ^7^* D 2 dA 


(35) 


Now we replace the square of the acceleration by its time average over a 
complete cycle of the motion: 


<a2> tt™’ <' 5 o“ 2co “' :>2) ti M 2 


fl 2 2 4 

V » 


(36) 


We note that the dissipation, being proportional to the entropy production, 
must be frame -independent, so that the value calculated in the instantaneous 
reference frame of the wall may be taken ever to the laboratory frame without 



change. Using (35) and (36), and integrating over the area of the entire cell, 
we get 


E 


.2 2 ,2 4 
D ra nolju 


LtiL 



R 


3 

r dr 


(37) 


where R is the radius of the cell; we then have 


v m 1 O 2 ” 2 " - 9jL 4 R 4 . -L. 

E 4 p_ V R x. 64rr np, 


M 2 a 2 4 
M 9^ 


H 


(38) 


where M is the total mass of the fluid in che cell. 

This relation connects the function (U)/x. which was calculated in section 
III and plocted in fig. 2 with the dissipation in the disc of fluid, the 
latter being proportional to the drive voltage, modulo a constant. The propor- 
tionality of (u)/k to this voltage enables us to see what the theory implies 
about the experiments. For purposes of comparison let us define x.^ as chat 

value of x. for which che derivative of (u)/x. with respect to x. has halt of 

2 2 

its large-R asymptotic value. Similarly we define T^ as that value of T for 

2 

which the derivative of the drive voltage with respect to T“ has half of its 

2 

limiting value. By expressing this T^ as a Knudsen number, we can determine 

whether the curvature in the experimental plot falls in the right temperature 

range to be the mean free path-finite size effect. 

In fact it does nor. We find that x.^ a 1 for all values of s. On the 

other hand the temperature T^ corresponds to values of k ^ 36 for p * 29 bar 

and x. — 7 for p ■ 0. Furthermore, we see from Table 2 that for a cell of the 

size used at Cornell, the curvature predicted by the theory should only be seen 

at pressures of a few bar, but the experimental curvature is seen at high pres- 

2 

sures as well. We conclude that the curvature in the drive vs. T plot is not 
due to the mean free path effect explored in this paper, and its explanation 
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raise be sought elsewhere, perhaps In the thermometry; the effect calculated here 
is being masked by this other, unknown effect. 

There are three courses of action that might be followed to actually see 
the mean free path effect, (i) Identify the cause of the curvature in the 
experimental plot, and remove it; this should unmask the upward curvature due 
to the mean free path effect, which should be noticeable at low pressures. 

(li) Remove the nuisance damping, so that the € -dependent vertical translation 
of the curves in fig. 2 can be determined experimentally; this part of the 
effect is significant even when the curvature is too small to be seen. (ill) 
Repeat the experiment in a thinner cell, so that a given value of K will corres- 
pond to a higher temperature. 

Finally, we note that in the superfluid phases the mean free path is 

4 

expected to become very large , so that almost all experimental cells will 

have dimensions of similar or smaLler magnitude. The calculation given here 

applies only to tne normal phase, but the method should be readily extensible 

14 

at least to the B-phase. We hope to treat this problem in a future publication. 
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Tab la I 


X. 6 


0.2 

.5869 

0.4 

.3401 

0.6 

.2266 

0.3 

.1609 

1.0 

.1176 

1.2 

.0868 

1.4 

.0638 

1.6 

.0469 

1.8 

.0330 

2.0 

.0225 


Deviation 5 of the function (U>/x. from its linear 
behavior at high x.. 


Table 2 


i , bar 

V mK 

X , um 

*, D*95um 

*, D-20um 

5 

0 

1.0 

68.0 

1.40 

0.295 

.4434 

3 

1.4 

29.2 

3.25 

0.684 

.1959 

6 

1.6 

19.5 

4.87 

1.03 

.1126 

12 

2.08 

9.36 

10.1 

2.13 

.0174 

18 

2.31 

6.26 

15.2 

3.20 


24 

2.45 

4.65 

20.4 

4.29 


30 

2.56 

3.64 

26.1 

5.49 



The mean free 

path and Knudsen 

number evaluated 

at the 



transition temperature for several pressures, and for two 
cell thickness. We have used data from Wheatley.^ In the 
last column, the deviation 6 is evaluated at the transition 
temperature, for a cell width of 20 jm. It is negligible 
at pressures above about 12 bar. 
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Figure Captions 


Figure l Sketch of Che idealized viscometer cell, showing an example of the 

velocity profiles found In the calculation. The curves represent 

u x (z) for x • 3 and •: • 0, 0.2, 0.4 and 0.6; each mark on the 

★ 

horizontal axis represents 5FT/m . 

Figure 2 The quantity (U)/x as a function of x for several values of •; . 

We have (U>/x ■ p ? (u)/FD where (u) Is the flow velocity averaged 
across the cell, F is the driving force per particle, D is the 
cell width, and p„ is the Fermi momentum. 

Figure 3 The ratio of the "effective viscosity" to the true shear viscosity 
of a bulk sample of the fluid, as a function of x, for several 


values of z 
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